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Abstract 

Tracy and Widom showed that fundamentany important kernels in random matrix 
theory arise from systems of differential equations with rational coefficients. More gener- 
ally, this paper considers symmetric Hamiltonian systems and determines the properties of 
kernels that arise from them. The inverse spectral problem for self-adjoint Hankel operators 
gives sufficient conditions for a self-adjoint operator to be the Hankel operator on L^(0, oo) 
from a linear system in continuous time; thus this paper expresses certain kernels as squares 
of Hankel operators. For suitable linear systems {—A, B, C) with one dimensional input and 
output spaces, there exists a Hankel operator F with kernel (p(x){s + ^) = Ce~^'^^~^^~^^^^B 
such that gx{z) = det(/-|- (2 — 1)FF^) is the generating function of a determinantal random 
point field on (0, 00). The inverse scattering transform for the Zakharov-Shabat system 
involves a Gelfand-Levitan integral equation such that the trace of the diagonal of the so- 
lution gives ^ log (7a; (z). Some determinantal point fields in random matrix theory satisfy 
similar results. 

Keywords: Determinantal point processes; random matrices; inverse scattering 
1. Introduction 

Traditionally, one begins random matrix theory by defining families of self-adjoint 
n X n matrices endowed with probability measures, known as ensembles, and then one 
determines the joint distribution of the random eigenvalues. By scaling the variables and 
letting n ^ 00, one obtains various kernels which refiect the properties of large random 
matrices. The kernels generate determinantal random point fields in Soshnikov's sense [16, 
20] . It turns out that many such kernels in random matrix theory have the form 



f{x)g{y) - f{y)g{x) 



x-y 

where / and g satisfy the system of differential equations 



{x,y > 0) 



d 

m[x) — 
ax 



'fix)' 




a{x) 


I3{x) ' 






_9{x) _ 




_-7(x) 


—a{x) 







(1.2) 



where m,,a,(3 and 7 are real polynomials. Tracy and Widom [17, 18, 19] began what 
amounts to a classification of kernels that arise from such differential equations, and their 
analysis revealed detailed results about the fundamental ensembles. 



This work was partially supported by EU Network Grant MRTN-CT-2004-511953 'Phenomena in High 
Dimensions'. 
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Of particular interest is the Airy kernel 

Aijx - X)Ai'{y - A) - Ai\x - X)Ai{y - A) 
Kx{x, y) = — — (1.3) 

on L^(0, oo) where Airy's function Ai satisfies Ai"{x) = a;Ai(a;). Some of the fundamental 
properties of this kernel involve the remarkable formula 

(>00 

Kx{x,y)= I Ai{x + u-X)Ai{u + y-X)du (1.4) 

which expresses the operator K as the square of the Hankel operator on L^(0, oo) that has 
kernel Ai{x + y — X). 

The differential equation (1.2) is an example of a symmetric Hamiltonian system, as 
we can define more generally. 

Definition (Symmetric Hamiltonian system). For an integer m > 1, let J be the matrix 



J = 



-/^ 

Im 



(1.5) 



which satisfies = —hm and = —J, and let E{x) and F{x) be (2m) x (2m) real 
symmetric matrices for each x > such that x E(x) and x i-^ F{x) are continuous. 
Then we consider the symmetric Hamiltonian system 

J^vE'A- (Ai?(a:)+F(a:))vE'A (1-6) 

where '^x{x) is a (2m) x 1 complex vector. In particular, when E{x) and F{x) have entries 
that are rational functions of x, we have a system considered by Tracy and Widom [19]. 
Given a solution G L°°{{0, oo); C^"^), we introduce the kernel 

x-y 

and we investigate the properties of Kg^x- 

More generally, we consider families of kernels Kt^xi^-, y) for A > 0, that satisfy some 
of the following properties as operators on H = L^(0,oo): 
(1°) the Lyapunov equation holds 

—Ks,x = -AKs,x - Ks,xA^ (1.8) 

as a sesquilinear form on D{A^) x D(A^), where e~^"^ is a bounded Cq semigroup on H 

and D[A^) is the domain of A^; 

(2°) < Ks^x < I for all s > sq for some sq < oo; 

(3°) s I— > Ks,A is decreasing and converges strongly to as s — > oo; 
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(4°) Ks^x is of trace class; 

ds 



(5°) the operator on H with kernel ^Kg^\ has finite rank. 



In fact, many of the properties of ensembles which arise in random matrix theory are 
essentially consequences of the properties (l'^) — (5*^), in a sense which we now make more 
precise. We recall from [16] the notion of a determinantal random point field. 

Definition (Configurations) . A configuration on R is an ordered list A = (Aj)°^_^ such 
that Xj < Xj+i for all j e Z; the configuration is locally finite if I'xi^L) = ^{j : Xj e L} is 
finite for all compact sets L. Let A be the space of all locally finite configurations on R. 
For each bounded and Borel set E, and A; = 0, 1, . . . , we let 

= {XeA:ux{E) = k} 

be the set of all locally finite configurations that have k terms in E; now let B be the 
cr-algebra generated by the . A random point field (P,A,S) on R is a probability 
measure P : S ^ [0, 1]. Wc let y{a,h) be the random variable that gives the number of 
points in (a, 6), so u^x, oo) = ^{j : Xj > x}. 



Definition (Correlation functions). Given nonnegative integers nj such that Xlj=i '"j 
and disjoint Borel sets Ej we consider A e A such that v\{Ej) > rij for all j. Then 

gives the number of ways of choosing rij points A^ from the v\{Ej) points of A that are in 
Ej for all j . The correlation function Rn : R"^ R+ for Pisa locally integrable function, 
which is symmetrical with respect to permutation of its variables, such that 

EiVp"j.,_1 u= I ••• / Rn{xi,...,Xn)dxi...dXn (1.10) 

1 fc 

for all disjoint Borel sets Ej {j = 1, . . . ,k). This is the expected number of configurations 
that have u\{Ej) > Uk for all j. 

Conversely, Soshnikov [16] observed that one can introduce a random point field from 
the determinants of kernels that satisfy minimal conditions. We state without proof the 
following general existence theorem for determinantal random point fields. 

Lemma 1.1. Suppose that K : R x R — > C is a continuous kernel such that 
(1°) the integral operator with kernel K{x, y) on I'^(R) satisfies < K < I; 
(4°) the kernel I[a,b]{^)^i^i y)^[a,b]{y) on L'^(R) is of trace class for all finite a, b. 
Then there exists a random point field such that the correlation functions satisfy 

Rn{xi,...,Xn) =det[K{xj,Xk)]2k=i (n = l,2, ...). (1.11) 
Further, the Rn (n = 1,2,...) uniquely determine P. 



Un = n 
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Definition (Determinantal point Geld). If the Rn have the form (1.11), then (P,A,B) is 
a determinantal random point field. 

In this paper we introduce natural examples of kernels K by means of linear systems, 
and recover properties (1°) — (5°) in a systematic manner. We summarise the construction 
in this introduction, and describe the details in section 2. Consider an operator A with 
domain D{A) in state space H such that the Co semigroup e~*^ is bounded, so ||e~*^ II ^ 
for some M < oo and alH > 0. Consider the linear system 

^ = -AX + BU (X(0) = 0) 
Y = CX (1.12) 

where B : C ^ D(A) is bounded and C : D(A) — C is admissible for A^ ^ soY ^ -^^(0, oo). 
We let (j){x) = Ce~^^B and 0(a;)(y) = (j){y + 2a;), then introduce the Hankel operators 

/■oo 

r<^(.)/(2/) = / (^i.){y + u)f{u)du (1.13) 

from a suitable domain in L^(0, oo) into L^(0,oo). We also consider the Gelfand-Levitan 
integral equation 

/>oo 

T{x,y)-^{x + y)- T{x, z)^{z + y) dz = {0 < x < y < oo) (1.14) 

J X 

where 5' and $ are both either (i) real scalars, (ii) 2x2 real diagonal matrices, or (iii) 
2x2 complex matrices. We state our main theorem as follows. 

Theorem 1.2. Suppose that the controUabihty Gramian 

poo 

L^= e-^^BB^e-*"^^ dt {x > 0) (1.15) 

J X 

is of trace class on H and of operator norm \\Lx\\ < 1; likewise suppose that the observ- 
ability Gramian 

nOO 

Q^= e-^^^C^Ce-^^dt {x > 0) (1.16) 



is of trace class on H and that \\Qx\\ < 1- 

(i) IfC = B^ and A = A^^ , then let g^iz) = det(/ + {z - l)T^^^^) and ^x) = (f){x). 

(ii) If (j){x) is real, then let gx{z) = det(/ -\- {z — l)r^^^^) and ^{x) = di?Lg[—(j){x), (t){x)]. 

(iii) Or let gx{z) = det(/ + {z - 1)T^^^^tI^^^) and 



(f){x) 
-(f>{x) 



{x > 0). 



Then in each case there exists a determinantal random point Geld on (0, oo) with generating 
function gx{z) = 'Ez''^^'°°^ such that 

d 

— log ^^(0) = trace T(x,x) {x > 0) (1.17) 

is given by the diagonal of the solution of the Gelfand-Levitan integral equation (1.14). 

The integral equation in case (i) is associated with the inverse scattering problem for 
the Schrodinger equation on the real line and in case (ii) by a pair of Schrodinger equations; 
whereas the integral equation in (iii) is associated with a Zakharov-Shabat system. 
The fundamental examples of determinantal random point fields in random matrix theory 
involve kernels associated with self-adjoint Hamiltonian systems of differential equations. 
In section 3 we introduce the notion of a symmetric Hamiltonian system with matrix 
potential, as considered previously by Atkinson and many others; see [5]. We consider 
spatial kernels Kx associated with symmetric Hamiltonian systems, and give a sufiicient 
condition for the kernel to factor as K\ = T^F^, where F;^ is a vectorial Hankel operator. 
As we show in section 4, this covers some fundamental examples of kernels that arise in 
random matrix theory, and we recover case (ii) of Theorem 1.2. A similar computation 
shows how (iii) arises. 

Schrodinger differential operators on L^(0, oo) with bounded potentials give rise to ker- 
nels in the spectral variables which satisfy (5°), as we discuss in section 5. The Korteweg-de 
Vries flow has a natural effect on the kernels. In section 6, we consider the Zakharov-Shabat 
system and establish case (iii) of Theorem 1.2; here the kernels behave naturally under the 
flow associated with the cubic nonlinear Schrodinger equation. Some of the calculations 
will be familiar to specialists in the theory of scattering from [1, 6, 23], but we include 
them here so that the paper is self-contained. 

2. Linear systems and their determinants 

Definition (Linear system). Let H he a separable complex Hilbert space, called the state 
space, and Hq a separable complex Hilbert space called the output space. Let e"*"^ be 
a Co semigroup on such that ||e~^"^|| < M for some M < oo and all s > 0, and let 
D{A) be the domain of the generator —A, which is a dense linear subspace of H and itself 

1/2 

a Hilbert space for the norm — (ll^lllf + ll^Cllif) • the language of linear 

systems from [14, 15], we consider the continuous time system 

^ = -AX + BU it> 0) 
Y = CX, X{0) = 0. (2.1) 

where B : Hq DiA) and C : D(A) — > Hq are bounded linear operators; this is known 
as a {-A,B,C) system. Let (j){x) = Ce-^^S, so G L~((0, oo); B(iyo)). The associated 
Hankel operator F<^ is the integral operator 

(>CXD 

T^fix)^ <P{x + y)f{y)dy, (2.2) 

^0 



5 



defined from some dense linear subspace of -£'^((0, oo); Hq) into -£'^((0, oo); Hq). 
We introduce the transfer function 



(l)iX) = CiXI + A)-^B, (2.3) 

which we recognise as the Laplace transform of (p{x) = Ce~^^B; the Fourier transform of 
(f) gives the scattering data. Suppose that U G L^{{0, oo); H) has Laplace transform U (A), 
and that : C-|_ — > B{Ho) is a bounded analytic function. Then Y e -^^((0, oo); Hq) has 
Laplace transform Y such that Y{X) = 0(A)C/(A). 

Definition (Admissible). We say that a bounded linear operator C : D{A) — > Hq is 
admissible for e~^'^ if Ce~^"^^ belongs to L^((0, oo); i^o) for all ^ e H, and there exists 
Kc{A) such that 

/•oo 

/ \\Ce-^''akd'<Kc{AmrH i^^H), (2.4) 
Jo 

equivalently, the operator : H ^ L^((0, oo); ii'o) is bounded where 0'''^ = Ce~*^^ 
and II ©11 = Kc{A). Examples in [9] show that the notion of admissibility is difficult to 
characterize simply. 

Definition (Schatten ideals). Let be the space of Hilbert-Schmidt operators, and 
be the space of trace class operators, with the usual norms. 

Proposition 2.1. Suppose that C is admissible for e~^^ and that B : Hq D{A^) has 
B^ admissible for e"*"^^ . Then the observability Gramian 

POO 

Qa:^ e-'^^C^Ce-'^^ds (x > 0) (2.5) 

J X 

and the controllability Gramian 

L^= I e-'^BB^e-'"^^ ds {x > 0) (2.6) 

J X 

define bounded linear operators Qx,Lx : H H by these strongly convergent integrals 
such that: 

(1°) the derivatives satisfy the Lyapunov equations 

^ = -A^Qx-QxA, ^ = -AL,-L,A^ (2.7) 

as sesquilinear forms on D{A) x D{A) and -D(A^) x D{A^) respectively; 
(2°) 0<Qx< Kc{Afl and < < Kst {A^fl for all x > 0; 
(3°) Qx and decrease strongly to zero as x increases to infinity. 

[4°) Suppose further that C{iyl + A)~^ is a Hilbert-Schmidt operator for all y ^ R and 
that ||C(%/ + A)~^\\'^2dy < oo. Then is trace class for each a; > 0, and 

tracego = — / \\C{ixI + A^^Wl^ dx. (2.8) 
27r 7_oo 
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(5°) Suppose that Hq has finite dimension m. Then rank^^ < m. 

Proof. (2°), (3°) The integrals converge by the definition of admissibihty, and the other 
statements are immediate consequences. 

(1°) For ^ G the e~*^^^ is diffcrentiable in H with derivative —e~^^^A^^. By 

the fundamental theorem of calculus, we have 

-AL^ - L^A^ = y ^ (e-"^SSte-^^^ ) ds = -e-^^SSte"^^^ , (2.9) 

as bilinear forms on D{A^) x D{A^), hence the result. 

(4°) Let {ej)°^i be an orthonormal basis for H. By Plancherel's formula in Hilbert 
space, we have 

pOO -1 /"OO 

\\Ce-y\j\\%Jy=—J_JC{ixI + A)-'ej\\%^dx (2.10) 
and summing this identity we deduce 

5^(Qoe,-,e,) = — / ^||C(za;/ + ^)-^e,||l,^dx (2.11) 



and hence 



1 

traceQo = WQTWI^ = 7^ / \\C{txI + Ay^f. dx, (2.12) 

27r 7-00 



so Qo is trace class. 

(5°) From the expression 



dx 



'C+Ce-^^ (2.13) 



it follows that the rank of is less than or equal to the rank of C, hence is less than or 
equal to m. 

□ 

Proposition 2.2 (Determinant of the observability Gramian). 

Suppose that (4°) holds, so that the observabihty operator : L^((0, 00); i^o) — ^ H is 
Hilbert-Schmidt, where 

poo 

ef= / e-'^'c^f{s)ds {feL\{0,ooy,Ho)). (2.14) 

(i) Then 

det(/ - XQ^) = det(/ - AP(^,oo)eteP(^,oo)) (A e C, x > 0). (2.15) 
deGnes an entire function that has all its zeros on the positive real axis. 
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(u) Suppose further that A — . Then the Hankel operator on L^((0, oo); i^o) with 
kernel ^{s + t) = Ce-(*+^)^Ct, has = G+G > and 

d 

— trace Qx = — trace c!>(2a;). 
ox 

(in) Suppose still further that Hq — C"^ where m < oo. Then the zeros of det(/ — XQx) 
have order less than or equal to m. 

Proof, (i) We have Q^i{t) = Ce'^^i and hence GGt = Qq- Further, since the operators 
G and G^ are Hilbert-Schmidt, we can rearrange terms in the determinant and obtain 

det(/ - AP(^,oo)0^0i'(x,oo)) = det(/ - AGP(^,oo)0^) = det(/ - AQ^). (2.16) 

The zeros of det(/ — AQa) are 1/Aj, where A^ are the positive eigenvalues of Q^- 

(ii) Now 

/•oo 

e^Qf{t) = Ce-'^ e-'^^C^f{s)ds, (2.17) 
Jo 

so G^G reduces to a Hankel operator when A = A'^ . Further, we have 

/•oo 

trace (5 a; = / trace e"*^C"' Ce"*^ <it 

J X 

/•oo 

= / trace Ce-'^^^CUt, (2.18) 

J X 

whence the result. 

(iii) The (block) Hankel operator with kernel 4)[x){s + t) — C'e^^+^+^^^'^C^^ is non 
negative and compact, and is unitarily equivalent to the some matrix [ajj^-k\Yk=i which is 
made up of m x m blocks. Hence its spectrum consists of together with a sequence of 
eigenvalues A^ of multiplicity less than or equal to m which decrease strictly to as j — > oo 
by [14, Theorem 2]. Hence the zeros of the function det{I — \P(^x,oo)^'^ ^P{x,oo)) have order 
less than or equal to m at the points 1/Aj. 

□ 

To express the hypotheses of Proposition 2.2(ii) in terms of spectra, we present the 
following result, which is known to specialists. 

Definition (Carleson measure). Let fi be a positive Radon measure on C4. — {z & C : 
> 0}. Then is a Carleson measure if there exists cq > such that 

IJ,{[0,x]x[y-x,y + x]) <cox {x>0,ye'R). (2.19) 

Proposition 2.3. Suppose that A is self-adjoint and has purely discrete spectrum (Kj), 
with Kj > listed according to multiplicity, and that {ej) is corresponding orthonormal 
basis of eigenvectors. Let (l){x) = Ce~^-^C^ . 

(i) Then T^j, is bounded if and only if \Cej\'^5^j is a Carleson measure. 

(ii) IfYl'jLi |C'ejp/Kj converges, then F,^ is of trace class. 
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Proof, (i) We use hats to denote Laplace transforms, and let H'^ be the usual Hardy space 
on C_|_ as in [10]. By the Paley-Wiener theorem, the Laplace transform gives a unitary 
map from L^(0, oo) onto H^. Then 

r-OO t'OO Q° 

{T^f,f)= / J] |Ce,f e-(^+*)'^V(^)M^^^^^^ 
Jo Jo 

oo 

= ^|Ce,f|/K)|^ (2.20) 

Hence F,^ is bounded if and only if there exists ci such that 

{J2\Cej\'S^^Afn<ci lim / \f{x + iy)\^dy; (2.21) 

which holds if and only if we have a Carleson measure; see [10]. 

(ii) Note that ■s/2k]/{z + kj) is a unit vector in i?^, so / i— >• 2Kjf{Kj)/{z + Kj) has 
rank one and norm one as an operator on H^; hence the result by convexity. 

□ 

Definition (Balanced system). If Qq = Kq, then the system is balanced. 
Remarks 

(i) The controllability operator S : L^{{0, oo); Hq) — > H 

/•oo 

S/= / e-'^Bf{t)dt (2.22) 
Jo 

satisfies an obvious analogue of Proposition 2.2. Note that F,^ = Q^E. 

(ii) One can interchange the controllability and observability operators by interchanging 

{-A, B, C) ^ (-A+, C^ 5t), which interchanges ^ fJ,. 

However, we will consider in section 5 some self-adjoint Hankel operators which do not 
arise from the special case of A = and B = . 

(iii) The set K of kernels that satisfy (2°), (3°), (4°) and (5°) is convex; further, for K ^ K. 
and U e B{H) such that \\U\\ < 1, we have U^KU e K. 

(iv) If (1°) holds with a finite-dimensional H, then (5°) holds. But (5°) is not implied by 
finite-dimensionality of Hq. 

(v) For X > 0, the shifted system {—A,e~^'^B,Ce~^^) has observability operator Q^, 
controllability operator and, with (j)(^x){t) = Ce^^^^+^^^S, the corresponding Hankel 
operator is F^^^^ . 

Proposition 2.4. (Determinants involving the Hankel operator) Suppose that the con- 
troUabihty operator and the observabihty operator 'E.^ for {—A, e~^'^B, e~^^C) are 
Hilbert-Schmidt. Then the operator : H H, defined by 

/"OO 

Rx^= / e-y^BCe-y^'^dy, (2.23) 

J X 
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is of trace class and satisGes 



det(/ - Ar<^(^) ) = det(/ - \R^). (2.24) 

Proof. By Proposition 2.2, the operator is trace class. By rearranging, we obtain 

det(/-Ar^(^,)=det(/-AetS,) 
= det(/-AS,et) 

= det(/- Ai?a;). (2.25) 

□ 

In section 3 we introduce some kernels that arise from Hankel operators as in Propo- 
sition 2.2. The kernels are defined with symmetric Hamiltonian systems, as we recall in 
section 3. 

3. Kernels arising from Hamiltonian systems of ordinary diflferential equations 

Let Z) be a domain that is symmetrical with respect to the real axis and contains 
(0, cxo). Wc later define kernels Kx that satisfy some of the following properties: 
(6°) Kx defines a bounded linear operator on H for all A e D; 
(7°) A 1-^ Kx is analytic on D; 
(8°) Kx = kI for all X e D; 

(9°) Kx is a Hilbert-Schmidt operator for all X E D; 

(10°) Kx is an integrable operator on L^(/; dx), for some interval /; so there exist locally 
bounded and measurable functions i/jk and such that 



m 



and J2T=i V'jl^; A) = 0. 

Definition (Hamiltonian system). Let E{x) and F{x) be (2m) x (2m) real symmetric 
matrices for each x > such that x i— > E{x) and x F[x) are continuous. Then we 
consider the symmetric Hamiltonian system 

J^*a(^) = {XE{x) + Fix))^xix) (3.1) 

where "^xix) is a (2m) x 1 complex vector. Suppose that for some A G C, the solution ^x 
of (1.4) belongs to L°°((0, oo); C"^). With the bilinear form {{zj), {wj)) = Xl^^i zjWj, let 

K4x,,)=S^HM^) = <£M£)l^^. (3.2) 

x — y X — y 

as in FHopital's rule, the diagonal of the kernel is taken to be 

Kx{x,x) = {^x{x),iXE{x) + F{x))^xix)). 
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Proposition 3.1. (i) Suppose that the solution belongs to L°°((0, oo); C^"^) for all 
X E D. Then K\ defines the kernel of a bounded linear operator K\ on I'^((0, oo); C) such 
that K{ = K-^, so (6"), (7°), (8°) and (10°) hold. 

(a) Suppose that E and F are bounded and that \E'a is a solution in L^{{0, oo); C^"^) of 
(3.2). Then Kx{x, y) defines a Hilbert-Schmidt kernel on i^^((0, cxo); dx), so (9°) also holds. 

Proof, (i) Indeed, the Hilbert transform H with kernel l/{n{x — y)) is bounded on L'^(R) 
and Kx is a composition of H and bounded multiphcation operators. The conditions 
(6°), (7°) and (8°) follow from basic facts about differential equations as in [8]. 

Observe that (J\E'a(x), "^fxix)) = 0, since we have the bilinear product; so the formula 
for Kx extends by continuity to a continuous function on (0, oo)^ and Kx is an integrable 
kernel as in (10°). 

(ii) There exist constants ci and C2 such that ||£?(a;)|| < ci and ||F(a;)|| < C2; hence 
the differential equation gives ||\&';)^(x)|| < (ci|A| +C2)||\I/A(a^)||- We deduce that belongs 
to L^{{0, 00); C^"^), and it is then easy to see that \I/a is bounded. A further application 
of the differential equation shows that ^'^ is also bounded. 

We split the Hilbert-Schmidt integral as 

Jo Jo J J\x-v\>i F ~ y\ 



+ 



x-y\<l 



j^,^,^^lM^lM)-d,dy. (3.3) 
x-y I 



The preceding estimates show that both of these integrals converge. 

□ 

Definition (Shift). For t > 0, let St : L2(0,oo) ^ L'^(0, 00) be the shift Stfix) = f{x-t), 
so that SlSt = /, and StSl = P{t,oo) is the orthogonal projection onto L^[t, 00) c i^^[0, 00). 
In the remainder of this section, we are concerned with the effect of the shift on solutions 
St : "^xix) ^ "^xix — t), and consequently on the kernels Kt^x = SjKxSt- 

Definition. For / an interval in R, let L : / x / ^ M^(C) be a continuous function. We 
write L >: if there exist continuous functions Ej : I ^ M^(C) for j = 1,2,... such that 

00 

L{x,y) = J2Ej{yyE^{x) {x,yel). (3.4) 



Lemma 3.2. Let v be probability measure v on I, and suppose that L y on I x I. 

(i) Then $ : L'^{iy; C"^) x L^(z^; C"^) — > C, is a positive sesquilinear form, where 

mv)- ff {L{x,y)ax),v{y))cMdxMdy) {tv^L\l;u;C^)). 
J Jlxl 

(ii) Suppose further that the defining sum (3.4) for L is finite. Then $ : L^(i/; C) 
L^(i/; C"^) has finite rank. 



11 



Proof, (i) The kernel is positive since 

oo „ „ 

= ^^'^k,i{^ I Ej{x)^k{x)i^{dx), / Ej{y)it{y)p{dy)^ > 0. 
k,e 3=1 k,e -'^ -'^ 

(ii) This is clear, since $ can be expressed a finite tensor. 

□ 

Proposition 3.3. Suppose that for some A > there exists a bounded and continuous 
solution '^x e L'^{{0-i oo); R^"^) to (3.2), where the coefficients E{x) and F{x) are bounded 
and satisfy 

m^m^Q, EM^m^o (x,,e(o,oo)). 0.5) 

X — y X — y 

(i) Then Kt,x = SlK^St satisGes (1°), (2°), (3°), (6°),(7°), (8°), (9°) and (10°). 

(ii) The kernel Kt^x is of trace class, as in (4°) and satisfies 

nOO 

tmceKt,x = {^x{x), {>^E{x) + F{x))^x{x)) dx. 

(Hi) If the sums involved in (3.4) for {E{x) - E{y))/{x - y) and {F{x) - F{y))/{x - y)) 
are finite, then (5°) also holds. 

Proof, (i) We observe that if L ^ 0, then {L{x, y)^, ^) gives the kernel of a positive definite 
operator. Kx{z + t,w + t) gives the kernel that represents SjKxSt, and hence satisfies the 
Lyapunov equation 

-Kt,x = -AKt,x - Kt,xA^ 

where —A — generates the semigroup 5'/. 
Prom the differential equation, we have 

d , ^ , / E(^x + t) - E(y + t) ^ , , ^ , 
g-^Kx{x + t,y + t)^- '-"^xix + t), ^x{y + t) 

F(x + t) - F(y + t) ^ , , ^ , ,\ 

^ ^ ^ ^-^xix + t), *A y + * ), 3.6 

X — y I 

and by the hypotheses on E and F we deduce that there exist Ej,Fj : (0, oo) M2m{C) 
such that 

Kx{x + t,y + t) = - Xj^i^ji^ + + t), Ej{y + t)^x{y + t)) 



dt 



Y,{Fj{x + t)^x{x + t), F,{y + t)^x{y + t)). (3.7) 
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The right-hand side gives the kernel of a negative definite operator, so 

^^{S}KxSJJ)<0 (3.8) 

for all / e L2(0,oo). 

By arguing as in Proposition 3.2, we see that Kx : L^(0, oo) L'^{0, oo) is compact. 
For / G iv^(0, oo), we observe that Stf weakly as t ^ oo and since Kx is compact 
K\Stf ^ in norm as t — oo; hence {S\ KxStf, f) 0. Consequently {Kt^xf, f) decreases 
to as t — > oo. Further, SjKxSt — in Hilbert-Schmidt norm as t — > oo, so there exists 
So such that ||5'/KA'S't || < 1 for all s > sq. 

(ii) We have proved that Kt^\ > 0, so the kernel is positive definite and continuous. 
By Mercer's trace formula, 

nOO pCO 

tTaceKt,\= / Kt,\{x,x) dx = / Ko,x{x,x)dx. (3.9) 
Jo Jt 

(iii) If the sum over j has finitely many terms, then the expression for ■^Kx{x+t, y+t) 
is a finite tensor and hence a finite-rank operator. 

□ 

We now relate the notion of positivity from the previous definition to matrix monotonicity 

in Loewner's sense. 

Definition (Matrix monotone) . Let / be an open real interval and let /'^ = R\/. Suppose 
that E : C\I^ ^ Mm{C) is an analytic function such that E{x) = E{x)'^ for all a; e / and 

{E{z) - E{z)^)/{2i) = QE{z) > {Qz> 0). (3.10) 

Then i? is a Loewner matrix function on /; equivalently, E is said to be matrix monotone. 

Theorem 3.4. Suppose that for some e > the functions z ^ E{z) and z i— > F{z) are 
matrix Loewner functions on (— e, oo). Suppose further that for > there exists a 
bounded and continuous solution G L^((0, oo); R^"^) to (3.1) such that "^fxi^) — > as 
X — > oo. 

(i) Then for all 3fJA > there exists 4> '■ (0, oo) — > Hq and a bounded Hankel operator 
T^^ : L2(0,oo) L'^{{0,oo);Ho) such that 

Kx = rlT^„ (3fJA>0) (3.11) 

and the family of kernels Kt^x — SlKxSt, for ?R.X > 0, satishes conditions (l°)-(4°), and 
(6°) -(8°), (10°). 

(ii) If Hq has finite dimension, then Kt,x aiso satisfies (5°). 

(iii) If Hq = C, then F<^^ is a scalar Hankel operator and Kx = F^^F,^^. 

(iv) If Hq = C and (f)x is real-valued, then Kx = F? . 
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Proof, (i) We need to obtain a suitable for the vectorial Hankel operator. Let D = 
C \ (— oo,£], and let Kx{z, w) be a kernel on D x D. Then K\{z + t,w + t) gives the kernel 
that represents SjKxSt, and Kx = F^^F^^ if 

(^)^^^i^A(t + z,t + w) = (M^), 4>xH)- (3.12) 

We have, from the differential equation, 
f d d , , / E(x) - Eiy) , , , , , A / Fix) - Fiy) , , , , , A 

(3.13) 

By the hypotheses on E and F, there exist self-adjoint matrices Ei,Fi > 0, self-adjoint 
matrices Eq and Fq, and positive matrix measures Qe and fii? such that 

= + + (3.14) 

and 

= Fix + Fo + / (— ^ - ^— (3.15) 

hence (3.13) equals 

^ E{x)-E{y) F{x)-F{y) 
X — y X — y 

=-,E.-,r ^ , -F.-r , (3.16) 

(w + a;)(w + 2/) (M + x)(w + y) 

By a straightforward Hilbert space construction similar to that in [3], we can introduce 
Hq and (j) e F^((0, oo); Hq) such that 

+ -^)[Kx{x,y)) ^ -{4>x{x),4>x{y))Ho- (3.17) 



Hence we have 
and 

and hence 



+ ^) (^^(^'^^ " rt^r,,(a;,y)) = 0, (3.18) 

Kx{x, y) - V\V^^ (x, y) ^ (x, y ^ oo) (3.19) 

/>oo 

KA(a;,y) = rJ,^r<^,(a;,y)= / (0A(s + x),0A(s + y))ds. (3.20) 

Now is bounded since Kx is bounded. Then for any Hankel operator S'^F,^^ = T^^S't. 
Hence, ^/fJ^^F^^^^ = V\^P^t^^)V^^ so that 

= SlKxSt = fJ.^P(,,^)F^, (A > 0). (3.21) 
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(ii) When Hq has finite dimension, the kernel {(j)\{x),(j)\{y)) has finite rank by Lemma 
3.2(ii). 

(iii) When Hq = C, the kernel of the Hankel operator is scalar- valued. 

(iv) In particular, the Hankel operator with (f)x '■ (0, oo) — > R is self-adjoint and 

□ 

Corollary 3.5. Suppose that K\ = F^^F^^, as in Theorem 3.4. Then under the Laplace 
transform C : L'^{{0, oo); C^"^) H'^{C+; C^"*), the nuUspace of K is unitarily equivalent 
to 6'_f/'^(C_|_; C^"^) for some bounded analytic function 9 : C+ M2m{C) that has unitary 

boundary values almost everywhere. 

Proof. The null space of K\ equals the null space of F,^^ and hence is a closed lin- 
ear subspace of L^((0, cxd); C^"^) which is invariant under the shift. Beurling's theorem 
characterizes the images of such subspaces under the Laplace transform; see [10]. 

□ 

Asymptotic forms of the differential equation as a;, A — > oo 

We now consider circumstances under which (3.2) does have a bounded or solution 
\&A- Suppose that E and F are is as in Theorem 3.4 and that Fi = in (3.14), so that F{x) 
is bounded. Then there are the following basic cases (i), (ii) and (iii) for the asymptotic 
form of (3.1) as A — > oo and x — > oo. 

(i) Suppose that Ei = 0. Then as a; — > oo we have E{x) — > Eq where 

r°° u 

Eo = Eo+ -—^VLE{du), (3.22) 

J e i- -r U 

and the asymptotic form of the differential equation is 

J^^x{x) = \Eq^x{x) (3.23) 

with solution 

^x{x) = exp(-AxJ^o)$A(0). (3.24) 

Now 3?(J£'o) = [JiEo\/2, so ^{JEq) is self-adjoint and has trace zero; hence 3?(J£'o) is 
either zero, or has both positive eigenspaces and negative eigenspaces. So in the following 
sub-cases, the solution of (3.23) has either: 
(i)(a) $A constant; 

(i)(b) $A oscillating boundedly as a; — > oo; or 

(i)(c) exponentially decaying solutions and exponentially growing solutions. 

In sub-cases (i)(b) and (i)(c) there exist bounded solutions $a to (3.23) such that 

{JMMm. (3.25) 

x-y 

gives a bounded linear operator on L^(0, oo). 
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(ii) Suppose that Ei is strictly positive definite. Then the asymptotic form of the 
differential equation is 

J^^A(a:) = XxE^^xix) (3.26) 

with solution 

^xix) = exp(-Aa;'Ei/2)$A(0), (3.27) 

and we have sub-cases (a), (b) and (c) analogous to those in (i) above. 

(iii) El of rank 1, . . . , 2m — 1. This case includes variants of Airy's equation. 

Examples 3.6. (i) Sonine considered the one-parameter families of functions Yi, that 
satisfy the system 



2u 

Yy-i + Yi^+i — — Y^^ 
z 



(3.28) 



as in [21, p 82]; the Bessel functions Yi, — Jj^ give solutions. One can transform the 
differential equation for Jjy into the system 



A. 
dx 



u 


= J 


' -l/x 


V 













-1 





u 




V 



(3.29) 



which has solution u{x) = \fx^v{^\fx). This system is matrix monotone when v = \. 

(ii) Theorem 3.4 applies to the Airy kernel (1.3), which describes the soft edge of 
certain matrix ensembles. Likewise, the Bessel kernel describes the hard edge; see [3, 4] 
for details. 

4. Determinantal random point fields 

In section 3 we showed how some important kernels factorize as K = r'l"r. Here we 
consider the properties of det(/— AF'^T); so first we introduce and solve the Gelfand-Levitan 
integral equation. 

Lemma 4.1. Suppose that —A : D{A) ^ H ^ H is a generator of a bounded Cq 
semigroup, B : Hq D{A), and C : D{A) Hq are linear operators, where Hq has finite 
dimension m, and let (j){x) — Ce~^^B. Suppose further that Kc{A), KBtiA"^) < 1. 
(i) Then the m x m matrix kernel 



Tx{x,y) = -XCe-'''^{l + XRa,) \-y^B (0 < a; < y, |A| < 1) 
gives the unique solution of the integral equation 

nOO 

Tx{x, y) + (l>{x + y) + X Tx{x, z)(j){z + y) dz = {0<x<y), 

J X 

and the kernel Tx{x, y) satisfies 

d'Tx d^Tx 



(4.1) 



(4.2) 



dx"^ dy^ 



q{x)Tx{x,y) = 



(4.3) 
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where q{x) = -2-^Tx{x,x). 

(a) Suppose further that m = 1, and that Qx ^nd are Hilbert-Schmidt. Then the 
determinant satisGes 

Tx{x,x) = ;^logdet(/ + Ar^(,^) {x > 0). (4.4) 



Proof, (i) First, we have ||-Ra;|| = ||Sa;Oj.|| < 1, so / + XRx is invertible and Tx{x,y) is 
well defined. One checks the identity by substituting the given expression for T\ into the 
integral equation. Further, ||r<^, , || = ||0j,S2;||, so / — AF ,t is invertible so hence solutions 

to the Gelfand-Levitan integral equation are unique. 

One can differentiate the integral equation and integrate by parts to obtain 

d^Tx d^Tx , , , , rfd^Tx d^Tx 



dx^ dy^ 
so by uniqueness 



+ Xq{x)cf>{x + y) + ^l [-9^- -9^) Hz + y)dz = 0, (4.5) 



(ii) When Hq — C the kernel takes values in C. Here — '^x^x trace class, and 
we can rearrange the traces and compute 

Tx{x,x) = -Atrace(Ce-^^(/ + Ai?a,)"^e-^^S) 
= -Atrace((/ + Ai?a;)"^e-^^5Ce-^^) 

= trace log (/ + XRx) 
ax 

= ^logdet(/ + AF^^^^), (4.7) 
where the last step follows from Proposition 2.2. 

□ 

Our first application is to the context of Theorem 1.2(i), where we consider determinantal 
random point fields associated with the observability Gramian. 

Theorem 4.2. Suppose that (4°) 6 : L^{{0,oo);C) H deGnes a Hilbert-Schmidt 
operator, and that (2°) the operator norm is ||0|| < 1. 

(i) Then there exists a determinantal random point held on (0, 00) such that ^{x, 00) is 
the number of points in {x, 00) and such that the generating function satisGes 

Ez''^'='°°^ =det{I+{z-l)Qx) {zeC,x>0). (4.8) 

(ii) Let F be the cumulative distribution function 
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Then 

F'{x)/F{x) = tj:ace{{A + A^)Q^{I-Q^)-^) (x > 0). (4.10) 

(iii) In particular, if A = A''', then det(/ + {z — 1)T^^^^ ) gives a generating function. 

(iv) When A = A''^ , the kernels 

Tx{x,y)^-XCe-'''^{I + XQ^)-^e-y^C^ (0 < x < y, |A| < 1) (4.11) 
and (j){x + y) = Ce~^^'^y^'^C^ satisfy the Gelfand-Levitan integral equation 

I'OO 

Tx{x,y) + X(l>{x + y) + X / Txix, z)(f)iz + y) dz = (4.12) 

J X 

and the diagonal satisfies 

Tx{x,x) = ^logdet(/ + Ar<^(^,). (4.13) 

Proof, (i) The kernel Ce~^^e~*"^^ of ©'''0 gives an integral operator on L^(0, oo) such 
that < ©'''0 < /; hence by Lemma 1.1 is associated with a determinantal random point 
field such that 

E^-(-'-) = det(/ + {z- i)eteP(,,^)) 

= dct(/ + (^-l)©P(,,^)©t) 

= det(/ + (^-l)g^), (4.14) 

so the determinant involving the observability Gramian gives rise to the determinantal 
random point field. 

(ii) We consider the probability that all of the random points lie in (0, x). The operator 
/ — Qx is invertible since ||(5a;||op < 1, and we have 

F(a;) = det(/ - Qa,) (a; > 0). (4.15) 

By a familiar formula for determinants, we have 

^logdet(/ - Qa;) = ^tracelog(/ - Qa;) 



/ d 
-trace(^(/ - Qx)~^-j^Qa 

trace((/ - Qx)~\A^Qx + Q^A)^ > 0, (4.16) 



where the last step follows from the Lyapunov equation (1°). Condition (3°) reassures 
us that F{x) is indeed an increasing function, and that the associated probability density 
function satisfies (4.10). 
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(iii) UA = A\ then Ce'^'+^^^C^ is the kernel of T^^^^ = Ot 6^;, so 

det(/ +{z- 1)Q^) = det(/ +{z- l)r^^J. (4.17) 

(iv) This is a special case of Lemma 4.1. 

□ 

Now we state the variant which arises in random matrix theory as in Theorem 3.4(iv) and 
Theorem 1.2(ii). 

Theorem 4.3. Suppose that A, B and C satisfy the hypotheses of Lemma 4.1, and that 

(j) = ,/)t. Let = ^(2x + y). 

(i) Then there exists a determinantal random point held on (0, oo) such that ^{x, oo) is 
tie number of points in {x,oo) such that the generating function satishes 

E^Ko^.oo) ^ ^^^^j ^ _ ;L)r2^^^ ) e C, o; > 0). (4.18) 

(a) Further, 

^logdet(/- A^r^^^^) =T_a(x,x)+Ta(x,x) (|A| < 1), (4.19) 

where Tx{x, y) — — ACe~*^(/ + XRx)~^e~y'^B satisGes a Gelfand-Levitan equation as in 
(4.2). 

Proof. First we check that = r|^^^ satisfies (2°), (4°) and (5°). We have Hr^^^J < 1, 

so < < /, and = Ol-Rx is of trace class. Hence by Soshnikov's theorem, we can 
form a determinantal random point field with generating function as above. 
To calculate the determinant, one can use the identity 

logdet(/ - X^K^) = logdet(/ - AF^^^,) + logdet(/ + AF^^^,), 

= log det(/ - XRx) + log det(/ + XR^). (4.20) 

(ii) The terms on the right-hand side satisfy 

^(logdet(/ - XRx) + logdet(/ + Ai?^)) = T_x{x, x) + Tx{x,x). 

by the Gelfand-Levitan equation. Indeed we can replace B in Lemma 4.1 by ±B to 
introduce iAF^^^^. 

□ 

We defer discussion of Theorem 1.2 (iii) until section 6. In section 5, we consider the 
determinant in Theorem 4.3 from the perspective of scattering theory. 

5. Scattering and inverse scattering 

The Gelfand-Levitan integral equation of Lemma 4.1 is closely connected to the 
Schrodinger equation, as we discuss in this section. Our aim is to identify a group of 
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bounded linear operators which acts naturally on the </> that appear in Theorem 4.3, and 
hence on the determinants. 

Given (j){x) = Ce~^^B as in Lemma 4.1, wc can solve the Gelfand-Levitan equation 

dx ' 



and recover q{x) — —24-T\{x,x). Further, given Tx as in Lemma 4.1, the function 



j>QO 

i;{x- k) = e^^^" + / e'y^T{x, y)dy (5.1) 

J X 

satisfies 

— -—^ip{x; k) + q{x)tl){x; k) — k^ip{x; k) (5-2) 

and 

il;{x;k)^e'^'' (x^oo). (5.3) 

This is a straightforward calculation, based upon (2.25). 

(i) Let be the discrete spectrum of — + g in L^(R), written \j = —Kj 
with Kj > so that each Xj = —Kj is associated with an eigenfunction i/j^x; Xj) that is 
asymptotic to e~'^^^ as x — > oo and Kn > ■ ■ ■ > ki > 0. We take c{—Kj) to be a constant 
associated with — 

(ii) The continuous spectrum is Eg = [0,oo), which has multiplicity two. For k E H 
and A = A;^ > 0, there exists a solutions ip{x; k) to (5.2) with asymptotic behaviour 



r.,M„-t*r'' !!! (5.4) 



e + 6(/c)e*^^ as X — >• oo; 
a{k)e~'^^^ as a; — > — oo. 

By [12], the reflection coefficient h belongs to Cq°, satisfies 6(0) = — 1 and b{—k) = b{k) 

and 

log(l-|6(fc)p) 



1 + k^ 



dk < oo. (5.5) 



(iii) By results from [12, 13], the transmission coefficient a extends to define the outer 
function 

1 p log(l-|&(t)p) 



on {k:Qk> 0} such that \a{k)\ = (1 - \b{k)\Y^^ for keK. 

The scattering map (j) associates, to the potential g, the function 

n- -1 /-oo 

(/.(x) = ci-K^fe-^i'' + — b{k)e''"' dk. (5.7) 

where the eigenvalues Xj = — ^j, the normalizing constants c(— k|) and the reflection 
coefficient 6(/t) are the scattering data. By (ii), (f){x) is real. 

The aim of inverse scattering is to recover q, up to translation from the scattering data. 
The following section uses computations which are extracted from [1,2], and originate 
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in other calculations of inverse problems as in [6,21]. The reflection coefficient b, the 
negative eigenvalues Xj and the normalising constants c{—Kj) determine q uniquely up to 
translation. 

Our general approach to the inverse spectral problem is to go 

(f)^ {-A,B,C)^T ^ q. (5.8) 

We now consider the first step in the process, namely realising {—A, B, C) from a given (f). 

Definition (Realisation). Given an bounded linear operator F, we wish to find a linear 
system {—A, B, C) such that the corresponding Hankel operator F,^ is unitarily equivalent 
to r. In particular, given scattering data (j){x), we wish to find a balanced linear system 
such that (^(x) = Ce~^^B for a; > 0. 

The following Lemma gives a characterization up to unitary equivalence of a special 
class of Hankel operators. 

Definition (Spectral multiplicity). For a self-adjoint and bounded linear operator A on 
H with spectrum S, let 

H= H{X)ii{dX) 
Js 

be the spectral resolution, where is a bounded positive Radon measure on S, such that 
Af{X) = A/(A). Now let 5{X) = dim.H{X) be the spectral multiplicity function for X E S. 

Lemma 5.1. Let T be a self-adjoint and bounded linear operator on H such that: 

(i) the nuUspacc of F is zero or infinite-dimensional; 

(a) F is not invertible; 

(ii) \6{X) - 5{-X)\ <1 for II almost all X. 

Then there exists a balanced linear system (—A, B, C) with Hq = C such that the Hankel 
operator on L^(0, oo) with kernel (j){x + y) = Ce^'^^+^^^S is unitarily equivalent to F. 

Proof. This is a special case of Theorem 1.1 on p. 257 of [14]. 

□ 

Proposition 5.2. Suppose that the hypotheses of Lemma 5.1 hold. 

(i) If the spectral density function b on the continuous spectrum is identically zero, then 
the scattering data can be realized by a linear system with finite dimensional H and Hq. 

(ii) Suppose that in the corresponding linear system A is a finite matrix such that all its 
eigenvalues Kj satisfy "Rkj > 0. Then the system is admissible. 

Proof, (i) By a theorem of Fuhrmann [15], one can choose ^ to be a finite-rank operator 
if and only if the transfer function is a rational function which is analytic on the closure 
of C+U{cx)}. 

(ii) The system is admissible since ||Ce~*^^|| < Me"***!!^!! for some M, k > and 

□ 

Remarks, (i) Not all self-adjoint Hankel operators satisfy the condition (ii) of Lemma 5.1. 
Consequently, there is a distinction between those self-adjoint Hankel operators that can 
be realised by linear systems in continuous time with one-dimensional input and output 
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spaces and the more general class that can be realised by linear systems in discrete time. 
In this paper wc concentrate on the continuous time case, while McCafferty has considered 
analogous results in discrete time, as in [11]. 

(ii) If 6 = in (5.9), then Propositions 5.2 and 2.3 apply to 4>. 

Definition (Evolution). For a system {—A, B, C), we refer to 4> as scattering data. Given 
a Co group E{t) on H and D{A), we can form the system {—A,B,CE{t)) and introduce 
E{t)(l){s) = CE{t)e~^"^B; thus the scattering data evolves with t. 

Article [3] highlighted the importance of groups that satisfy the Weyl relations; here 
we show that these are associated with evolutions that do not change the determinants in 
Theorem 1.2. 

Proposition 5.3. Suppose that D is skew sclf-adjoint and generates a Cq group of unitary 
operators that satisfies e^^^e~^^ — e-*"s*e^*"^e~^'^ for aii s € R, t > and some a € R. 
Then {—A, e^^B, Ce~^^) has observability Gramian Q^'^\ operator R^^^ and controllability 
Gramian L(«) such that (i) det(/ + AQ(«)), (ii) det(/ + Ai?(«)) and (Hi) det(/ + AQ(^)l(^)) 
are independent of s. 

Proof, (i) Wc have e'^"^^ e'^C^Ce-'^e-^^ = e'^e'^^^ C^Ce'^^e-'^ , so 
q{s) _ e^DQ^^-sD g^j^^ unitary equivalence det(/ + AQ^"'-') = det(/ + XQo). 

(ii) We have R^'^ = e-'^Roc'^, so det(/ + XR^'^) = det(/ + XRq). 

(iii) Likewise, L^^^ = c^^Lqc"'^^, and there is a unitary equivalence Q^^^L^^^ = 
e^^QoLoe-^^ leading to det(/ + AQ(*)l(*)) = det(/ + AQoLq). 

□ 

Lemma 5.4. Suppose that '^t{x; k) gives a differentiable family of vectors in C^"^ such 
that 

^M^; k) = Zt{x; k)-^t{x; k), (5.9) 



dt 
where 

at{x;k) /3t{x]k) 



Zt{x;k) 



--ft{x;k) -atix-k) 



(5.10) 



and at{x; k), Pt{x: k) and '^t{x; k) are symmetric m x m matrices. 

(i) Then with the bilinear form on C'^^ , the family of kernels 

Kt,x{i^,k) = {k,KeR,k^K) 

K — K 

satisfies (10") and 

k) = {j^^liBA^^yU^; k)) (5.11) 

(ii) If the at{x; k),j3t{x; k) and ^t{x', k) are rational functions of k, then -§iKt,x is of finite 
rank. 

Proof, (i) This follows by direct calculation, where the the effect of ^ is to replace J by 
JZt{x; k) + Zt{y; k)^J in the kernel. Then one uses the identities 

JZt{x; k) + Zt{x; kfj = J{Zt{x; k) - Zt{x; k)), (5.12) 

22 



which follow from the special form of the matrices. 

(ii) Given the formula (5.21), one can use the partial fraction decomposition of the 
entries to express -^Kt^x(K, k) as a sum of products of functions in the variable k or k. 

□ 

In accordance with the approach of [7], we are particularly interested in the case 
where k ^ Zt{x; k) is a polynomial such that the leading coefficient has trace zero. For 
Schrodinger's equation, we can introduce such families of matrices associated with the 
KdV flow. Let be the space of functions / : R ^ C that are infinitely differentiable 
and such that |a;p — as a; — cxd for j, £ = 0, 1, . . .. 

Suppose that q satisfies that q = v' + . Given a i/^ that satisfies —'4>" + qi/j = k'^i/jj 
we have a solution of the symmetric Hamiltonian system 

d 

dx 

Now let V evolve according to the modified Korteweg-de Vries equation 

dt dx^ dx ' 

and introduce functions of by 

« = - (l/2y , /9=(-l/2)(^, + i;2), ^={l/2){v,-v% 5 = v. (5.15) 







V ik 




i; 


p 




ik —V 




P 



(5.13) 



(5.14) 



Lemma 5.5. The matrices 

Vt{x;z) = 
give a consistent system 



V z 
Z —V 



and Zt{x; z) = 



^<f = Vtix;z)<i^, 
i^^Zt{x;z)^. 



a + bz^ I3z + z"^ 
7^ + z^ —a — 8z'^ 



(5.16) 



(5.17) 



Proof. As in [7] , it follows by direct computation that 

dVt{x;z) dZt{x;z) 



dt 



dx 



+ [Vt{x;z),Zt{x;z)]=0, 



(5.18) 



r^' and the system is consistent. The key idea is that one can equate 



dxdt ~ dtdx 

coefficients of the ascending powers of z, then one can eliminate the functions a,(3,'y and 
5 by simple calculus. 

□ 

Let ^t{x; k) be the solution of (5.17) that corresponds to z — ik where k belongs to 
R and k'^ to the continuous spectrum (0, oo). With the bilinear form (. , . ) on C^, let 



Kt,x{i^, k) = 



{J^t{x;K),^t{x; k)) 
i{K — k) 



(5.19) 
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where the numerator vanishes on k = so ^ is an integrable operator. This family of 
operators undergoes a natural evolution under the KdV flow, as follows. 

Theorem 5.6. Suppose that ^t(a:; k) give a locally bounded family of solutions which is 

differentiable in (t, x, k) and subject to \I/t(0; k) = for some G C^. 

(i) If v{x) — 0, then Kt^x is a multiple of the sine kernel [3, (1-2)]. 

(a) The kernels Kt^^ satisfy (10°) and (5°), so -S-Kt^x ^nd ■§TKt,x of finite rank. 



(Hi) The function u = ^ + v'^ satisfies the KdV equation, and as q{x) evolves to u(x, t) 

the scattering data for u undergoes a linear evolution (f) i-^ E{t)(f). 

(iv) Let {b, c{—K'j), Kj) be the scattering data for q{x), and suppose that b{k),b'{k) and 
k'^b{k) belong to -L^(R; dk). Then rE(t)4> gives a Hilbert-Schmidt operator for all t. 

Proof, (i) This is an elementary computation. 

(ii) Using Lemma 5.4, we calculate the derivatives, and find 



t,x 



dt 



-J + {k'^ + kK + k'^) 
iS{k + k) 



id{k + K.) 
+ kK + 



(5.20) 



which gives a kernel of finite rank, and likewise 



dx 



1 
-1 



*t(a;;K),*t(a;,A;) 



(5.21) 



which also gives a kernel of finite rank on L?{—oo, oo). 

(iii) By Miura's transformation, the function u = ^ -\- v'^ satisfies the KdV equation 



dt dx^ 



„ du 
ox 



(5.22) 



see [6, p. 65]. The evolution of the potentials u{x,0) h- > u{x,t) under the KdV fiow gives 
rise to a linear evolution on the scattering data. Now let ^t{x] k) be a continuous and 
uniformly bounded family of solutions of the system 



■£^t{x]k) = Ut{x]k)-<^t{x;k) 
d 



jr^t{x;k) ^Wt{x;k)^^t{x;k) 



(5.23) 



where 

Ut{x;k) 



1 
u-k'^ 



Wt{x;k) 



-1 



4z/c3 - fM 



2u + 4/c2 



(5.24) 

Then by considering the shape of the matrices in (5.24), we obtain the asymptotic forms 
of the solutions 



^'t(a;; k) x —ia{k)e 



-ikx 



"iftix; k) X —ie 



-ikx 



+ ib{k) 



I 
k 

—I 
k 



{x —>■ — oo). 



ikx-2ik^t 



{x oo); 



(5.25) 
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hence a(k) i— > a(k) and b(k) i— > b(k)e ^^'^^^ under the flow. By [6, p. 75], there is a group 
of hnear operators E{t) on the Hilbert space C" ® L'^(R) defined by 



E{t)(f>{x) = J2 c(-«f )^e-'^^^-2«?* + ^ h{k)e'^''- 



dk 



(5.26) 



such that u{x,y) corresponds to E{t)(p, and = max{e~^*'*", 1}. 

By applying Fourier inversion to the definition of the Airy function, we can express the 

integral over the continuous spectrum as 



1 r°° 



,ikx—2ik'^t 



-1 



(6t)i/3 J, 



Ai 



x-y 



{-6t) 



(iv) The Hankel operator with kernel J2^=i c(-K^)^e~''''^*+^^~^'^^ is clearly of trace 
class, so we need to consider the Hankel arising from b. To show that x{E(t)(/}{x))'^dx < 

oo, it suffices by Plancherel's theorem to show that b{k) and ^(6(A;)e~^*'^^*) belong to 
L^(R; dk). This follows directly from the hypotheses. 

□ 

6 Determinantal random point fields associated with the Zakharov— Shabat 
system 

In this final section we prove the remaining case (iii) of Theorem 1.2, and then we address 
the corresponding scattering theory. 

Consider the matricial Gelfand-Levitan integral equation 



(>00 

T{x, y) + \^{x + y) + A / T{x, z)^{z + y)dz = (0 < a; < y), 

J X 



(6.1) 



where, suppressing the dependence of T upon A, we write 



T{x,y) = 



U{x, y) V{x, y) 
-V{x,y) U{x,y) 



$(x) 



(f){x) 
-0(a;) 



(6.2) 



Theorem 6.1. Suppose that the system {—A,B,C) has Hq = C and with (f)(^^^{y) = 
^^-(2x+y)A^ satisfies, as in Lemma 4.1: 
(2°) lie^ll < 1 and < 1, and 
(4°) and Sa; are Hilbert— Schmidt. 

i) Then there exists a determinantal random point Geld on (0, oo) such that ^{x, oo) is the 
number of points in (x, oo) and such that the generating function satisfies 



g,{z) = E^-(-'-) = det(/ + {z- l)r^(.,r;^^^) 



(6.3) 



(a) Further log gx{z) = 2U{x,x), where U is given by the diagonal of the solution of 
the Gelfand-Levitan equation (6.1). 
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Proof, (i) We have < Qx ^ I and < < I. Let = P(^x,oo)'^^ LxQP{x,oo)i which 
defines a trace-class kernel on L'^{x,oo) and satisfies < Kj; < I. Then by Lemma 1.1, 
there exists a determinantal random point field on (x, oo) with generating function 

det(/ +{z- l)Kx) = det{I +{z- l)QP^x,oo)Q^ L,) 

= detiI+{z-l)QxLx), (6.4) 

and we continue to rearrange this, obtaining 

det(/ + {z- l)Kx) = det(/ + {z - l)P(,,^)0tSP(,,^)St0) 

= det(/ + {z - l)P(3.^oo)r<^-P(x,oo)r0-P(a;,oo)) 

^det(/+(^-l)r<^(^,rt^^^). (6.5) 

(ii) This identity is proved in the following two Lemmas. 

Lemma 6.2. Suppose that Hq — C and let 0(x) = Ce~'""^B and Gx — I + X^QxLx- Then 
the Gelfand-Levitan integral equation (6.1) reduces to 

POO foo 

V(x,y) + X^{x + y) + / V(x, s)(l){s + z)^(y + z)dsdz = (0 < a; < y), (6.6) 

Jx J X 

which has solution 

V{x, y) = -AS+e-^'"G-^e-^'^C^ (0 < x < y), (6.7) 
U{x,y) = X V{x,zmz + y)dz. (6.8) 

J X 

Proof. Once we have V, we can introduce U via (6.7), and the resulting matrix T satisfies 
the Gelfand-Levitan integral equation. To verify the equation for T, we first check that Gx 
is invertible when 3?A^ > — 1. The operators Qx and Lr^ are Hilbert-Schmidt and positive, 
so the operator QxLx is trace class, and hence the determinant satisfies 

det Gx = det(/ + A^Q^^L^Qy^) > o (6.9) 

since 3?(/ + X^QI^^ LxQI^^) > (1 - ^X^)I. 

One can postulate a solution of the form V{x, y) — X{x)^e~^^'^C^ , for some function 
X : (0, oo) — > H and by substituting this into the integral equation, one finds that X 
should satisfy 

X{x)^e-^^yC^XB^e-^^^''+y^C^ 

POO poo 

+AM / X{x)^e-'^^'C^Ce-^^'+^^BBh-'^^^'+y^CUsdz = 0, (6.10) 

Jx Jx 

so we want 

X{x)\l + X^QxLx) + A^te-^^^ = 0, (6.11) 



26 



and we can make this choice since is invertible. 



Lemma 6.3 The diagonal of the solution satisGes 



□ 



Proof. Prom (6.8), we have 

/'CO 

U{x, y) = -A^ / Ste-^'^G-^e-^'^CtCe-^^^+^^Sd^ 

J X 

= -X^B^e-^^'^G-^Qa^e-^yB. (6.13) 

Hence we can write 

= -A2trace(G-^ga.^L^). (6.14) 

We temporarily assume that A is real to derive certain identities, and then use analytic 
continuation to obtain them in general. Using Proposition 2.6, and rearranging various 
traces, we can derive the expressions 

Xhrace(G-^Qa,^^ = trace (^(GJ,)"^ A -A+ {Ga:)~^A^ - A^^ (6.15) 

and likewise 

AWe^G-i^L,) = trace ((G^)-^yl-yl+(G^)-^ylt_^t^^ (6J6) 
and since ^ = A2(^La; + Qc^^), we deduce that 

U{x,x) = ltrace{G-'^) 

1 , 1 ^ 
= o 3- log det Gx 

2 ax 

= li^^^^'^''^^ + ^"^^i^Aj- (6-17) 
This concludes the proof of the Lemma, hence of Theorem 6.1(ii) and Theorem 1.2(iii). 

□ 

We let q e Go°(R; C) and consider the Zakharov-Shabat system 



—ik q{x) 
—q{x) ik 



*(a;;A;) (6.18) 
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with '^{x; k) a complex 2x2 matrix. We observe that this matrix is skew-symmetric with 

zero trace, so the norm of any solution is invariant under the evolution, as is the Wronskian 
of any pair of solutions; hence the fundamental solution matrix of this system belongs to 
SU{2). We introduce the solutions '^+{x; k), "^-{x; k) e SU{2) such that 



^+{x;k) 



,—ikx 





,ikx 



[x — >• oo), 



(6.19) 



^-{x\k) 



-ikx 







-,ikx 



(x — > — oo); 



(6.20) 



then we introduce the scattering matrix S{k) e SU (2) such that ^'_(a;; k) — '^+{x; k)S{k) 
and we write 

"aik) P{k) ■ 
p{k) -a{k) 



S{k) 



(6.21) 



Now suppose that a and P are analytic on the upper half-plane, and that a has zeros at 
Kj. As in [6], we introduce the scattering data 



f^i^) Jkx 



{k) 



dk. 



(6.22) 



The sum contributes a function that decays exponentially as x — >^ cxo. 

Proposition 6.4. Let (—A, S, C) realise the scattering data (j) of the ZS system, suppose 
that the Gramians and are Hilbert-Schmidt. Then the potential satisfies 



1 cf 
2dx^ 



(6.23) 



Lemma 6.5. Let T be as in Lemma 6.2 and (6.2), and let 



*(a;; k) 



/"OO 

+ / T{x,y) 

J x 



^^iky 



dy. 



(6.24) 



Then 



^2 



dx"^ 



^{x; k) + W{x)'^{x; k) = k'^'^{x; k) 



(6.25) 



where W{x) = -2-^T{x,x). 

Proof. One can follow the proof of Lemma 4.1 and deduce that 

^^(^' y) - ^^(^' y) = W{x)Tix, y). 
Then one can verify the differential equation for '^{x; k) by direct calculation. 
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(6.26) 



From the original differential equation (6.18) we have 



*(a;; k) + 



-m Q 

-q' -\q\ 



'^{x;k) = A;^*(a;; k), 



so by equating the matrix potential with —2-^T(x,x), we obtain 



I |2 / 

— \q\ q 

-q' -IqI 



-2- 



dx 



U{x,x) V{x,x) 
—V{x,x) U{x,x) 



(6.27) 



(6.28) 



□ 



Finally, we consider how the potential q{x) evolves to u{x,t) under the nonlinear 
Schrodinger equation. Suppose that 



Wt{x;0 = 



-u iC 



' dx 



(6.29) 



Proposition 6.6. Suppose that u satisfies the nonhnear Schrodinger equation 

du d'^u ^1 ,2 

(i) Then the pair of differential equations 

£^ = Wt{x;0'^ 
^* = Z,(x;C)* 



(6.30) 



(6.31) 



gives a consistent system. 

(a) Let ^t{x, C) be family of solutions to (6.31 ) with initial value '^t{0, C) = Then the 
kernels 

Kt,x{K,k) = (6.32) 

satisfy (10°) and (5°); so -^Kt^x and ■§iKt,x have finite rank. 

(Hi) The evolution of the potentials under the NLSE gives rise to a linear evolution on the 
scattering data. 

Proof. This is similar to that of Theorem 5.6. 

□ 
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